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ABSTRACT 

6o ; 

^\ . Observations combined with theoretical modeling of the light curve of the recently discovered 

\ supernova 1997D in NGC 1536 suggest that it might host a black hole formed in the aftermath 

^ ' of the explosion. We consider some observable consequences of a black hole in SN 1997D and 

I ^ . estimate the late-time accretion luminosity of the material which falls back onto the hole. We 

\ find that this luminosity, decaying with a characteristic power-law dependence on time, may 

' emerge above the emission of the envelope in just a few years. Its detection would thus provide 

unmistakable evidence for the presence of a black hole. 
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(<~^ ' Subject headings: Accretion — black holes — hydrodynamics — radiative transfer — relativity 

rn ■ — supernovae: individual (SN 1987A, SN 1997D) 
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0^ ' 1. Introduction 

■ 

The nature of the compact remnant formed in the aftermath of a supernova explosion depends on the 
Q . mass of the iron core of the progenitor star, on the amount of material that falls back onto the core and 

y-tj \ on the maximum allowed neutron star mass (Mmax)- The core is usually in the range 1.4-2.0 Mq, while 

c/3 ■ the actual value of A/max is not conclusively established yet. According to recent results on the properties 

. . , of dense nucleon matter, the maximum (nonrotating) neutron star mass falls in the interval 1.8-2.2 Mq 

' (Akmal, Pandharipande & Ravenhall 1998), although smaller values are not yet completely ruled out. 

■ Consequently, the amount of matter fallback plays a crucial role in determining the fate of the collapsed 

^ \ remnant. An investigation of the fallback process for stars with different masses and different explosion 

■ - - ' energies has been performed by Woosley & Weaver (1995). They find that stars with main sequence mass 

larger than ~ 25 Mq have rather massive cores (~ 2Mo) and undergo accretion of a significant amount of 
matter, so that the remnant mass is larger than 3Mq, probably leading to the formation of a black hole 
(BH), after the launch of a successful shock. Fallback of a smaller amount of matter, ^ 0.1 — 0.3 M©, may 
occur in explosions involving a progenitor of 19-25 Mq, although the actual value depends sensitively on 
the energy of the explosion. Even in this case, a low mass (~ I.^Mq) BH may form if the equation for 
nuclear matter is not too stiff (Brown & Bethe 1994). The fact that a bright [L ^ 10'^^ erg s~^) accreting 
neutron star was not detected in SN 1987A in the early years following the explosion led to the hypothesis 
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that a light underluminous BH may have formed (Chevalier 1989; Brown & Bcthe 1994; Woosley & Timmes 
1996). Woosely & Timmes (1996) first suggested that a signature of BH formation in the light curve 
would be a bright plateau that falls steeply to very low or zero luminosity because of the low radioactive 
emission from heavy nuclei present in the ejecta. Recently, Zampieri, Colpi, Shapiro & Wasserman 
(1998a; hereafter ZCSW) pointed out that a BH may leave a distinguishable imprint in the light curve, 
which could be revealed when the main sources of radioactive decay fade away (see also Zampieri, Colpi, 
Shapiro & Wasserman 1998b). Using a spherically symmetric, general-relativistic, radiation-hydrodynamic 
Lagrangian code, ZCSW showed that fallback may lead to an accretion luminosity decaying secularly with 
time as L oc For SN 1987A, the accretion luminosity turns out to be too dim to be detectable, at 

present, above the radioactive decay emission of ^"^Ti. 

Recently, a new supernova, SN 1997D, was discovered in NGC 1536 (De Mello & Benetti 1997). 
The very low absolute luminosity and the low inferred heavy element abundance in the ejecta make this 
supernova a promising candidate for the detection of a compact remnant. Following the suggestion by 
Turatto et al. (1998), we consider the consequences of a BH forming in SN 1997D. We give an estimate of 
the accretion luminosity and the time of its emergence above the radioactive decay emission of the envelope. 



SN 1997D was serendipitously discovered on January 14, 1997 by De Mello & Benetti (1997) during 
an observation of the parent galaxy NGC 1536. Although the discovery of this supernova was made 
after maximum phase, a number of closely spaced observations of NGC 1536 set an upper limit to the 
optical luminosity ~ 2 magnitudes lower than that of a typical Type H supernova. This result makes SN 
1997D possibly the least luminous Type II supernova ever observed. A comparison of theoretical models 
with observations seems to indicate that the epoch of the explosion probably occurred 50 days before 
the supernova was discovered (Turatto et al. 1998). During all the evolutionary stages, the bolometric 
luminosity of SN 1997D is orders of magnitude lower than that of the "prototype" SN 1987A (see Figures 
|l| and 1^). At around 70-80 days after the explosion, the emission starts to be powered by 7-ray energy 
deposited into the ejecta and, as shown also by SN 1987A, the decline of the tail of the light curve follows 
closely the decay rate of ^^Co, the decay product of ^^Ni. If all the 7 rays are thermalized, the light curve 
after few weeks is expected to follow the relation 



where = 1.4 x lO^'^ erg s^^, tq — 111 days is the decay time of ^^Co and Mco is the total mass of ^^Co 
present in the ejecta. Assuming for SN 1997D the same 7-ray deposition as in SN 1987A and using the 
ratio of the luminosities in the tails of the light curves, Turatto et al. (1998) estimate the mass of ^^Ni 
ejected by SN 1997D to be about 40 times less than that in SN 1987A, i.e. 0.001-0.004 Mq. Another 
peculiar property of SN 1997D is the low expansion velocity of the ejecta deduced from measurements of 
the shift of the absorption lines in the spectrum. It turns out to be ~ 1200 km s"-'^. 

Using a spherical, Newtonian, Lagrangian hydrodynamical code to simulate the explosion, Turatto 
et al. (1998) were able to produce a satisfactory fit to the light curve of SN 1997D. To fit the data, a 
progenitor model with a radius r = 2 x 10^^ cm and a total mass M = 2QMq was adopted. The explosion 
energy was taken to be £^ = 4 x 10*^° ergs and the low luminosity tail from radioactive decay is reproduced 
by the ejection of i\/Ni = O.OO2M0 of ^^Ni. This figure is in agreement with that estimated by rescaling the 



2. The peculiar Type II supernova 1997D 
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value of SN 1987A from the ratio of the luminosities of the tails. The small values of r and E/M account 
for both the observed low peak luminosity and expansion velocity. 

The isotope ^^Ni is produced in the deep layers of the star, at the outer edge of the "iron core". In 
order to eject only a small amount of ^^Ni, a considerable fraction of the innermost stellar material has to 
fall back onto the central compact core where it is photodisintegrated before emitting 7-rays. The usual 
way to account for this effect in hydrodynamic simulations of the explosion is to impose a so-called "mass 
cut" in the post-explosion structure of the star. In the 26M0 model of SN 1997D of Turatto et al. (1998) 
the mass cut has been placed at 1.8 Mq. This guarantees the ejection of a mass of Ni as low as 0.002 
Mq. The location of the mass cut implies that the mass of the central compact remnant is ^ l.SM©. Then, 
as suggested by Turatto et al. (1998), a BH could have formed and the low luminosity tail of SN 1997D 
might be a signature of the formation of such an object. 

In the following we will need an estimate of the abundance of '"'^Co and "^^xi in the ejecta of SN 
1997D. As an upper limit to the abundance of these two isotopes, we take the value inferred for SN 1987A 
(~ 5 X 1O~^M0) and reduce it by a factor of 40, similarly to what estimated for the abundance of ^^Ni. 



In a recent paper we have presented the first self-consistent, fully relativistic investigation of supernova 
fallback in presence of a BH (ZCSW). We computed the light curve from a hot expanding hydrogen "cloud" 
in which the emitted energy comes from the release of the heat residing in the original gas or generated 
by compression in the course of accretion onto a central BH. Radioactive energy sources were not included 
in these calculations. At the onset of evolution the cloud is optically thick and has homogeneous density 
Pq] it is set into homologous expansion with a velocity profile u = r/tQ, with Iq the expansion timescale. 
The maximum velocity of the ejecta is Vq = u{rout) = I'out/to, where Tout is the outer radial boundary. 
During evolution, the unbound shells propagate outwards while those with an outward velocity smaller 
than the escape speed fall back after turnaround, causing accretion on a secular timescale. Initially, the 
cloud is radiation-dominated and photons diffuse outwards from the expanding stellar envelope. Hydrogen 
then recombines and the "cloud" becomes optically thin. Depending on the initial parameters, a significant 
amount of internal energy stored up into the gas is liberated, giving rise to a bump in the light curve. The 
density and velocity profiles display a self-similar character and are unaffected by the propagation of the 
recombination front. After the emission maximum, the gas internal energy is nearly exhausted and the 
luminosity falls off abruptly. 

Soon after the recombination wave has propagated inward to the inner accreting region, the front stalls 
because the compressional heating due to accretion surpasses radiative cooling. From this moment on, the 
light curve is entirely powered by accretion: the central BH "becomes visible" . The numerical calculations 
show that, at this stage, the accretion rate M in the inner regions is nearly independent of radius. The 
evolution in the luminosity-accretion rate plane closely follows the analytic curve L = L{M) derived by 
Blondin (1986) for stationary optically thick hypercritical spherical accretion onto BHs: 



3. Late— time accretion luminosity 




(2) 
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Fig. 1. — Bolometric light curves for SN 1987A (from ZCSW). The luminosity L is plotted as a function of 
time t. The triangles show the observed bolometric luminosity of SN 1987A. The dashed lines represent the 
expected contribution from the decay of radioactive elements (0.07 Af© of '^^Co and 5 x lO^^M© of ^^Co and 
44rpj-j rpj^g dotted line denotes the calculated bolometric luminosity emitted via accretion onto a putative 
BH in SN 1987A. The arrow marks the time of the BH emergence. 
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Fig. 2. — Same as in Figure 1 for SN 1997D. The data for SN 1997D circles are taken from Turatto et al. 
(1998). The dashed lines represent the expected contribution from the decay of radioactive elements (0.002 
Mq of ^^Co and 1.25 x 10~^A'fQ of ^^Co and ''^Ti). Curves a, b, c and d denote the calculated accretion 
luminosity at late times for the following values of the initial density and expansion timescale: a,) to — 2 
hr, po — 0.1 g cm~^, b) to = 2 hr, po = 2 x 10^^ g cm^^, c) to = 1 hr, po = 0.1 g cm"^, d) to — 2 hr, 
po = 2 X 10~^ g cm~^. The continuous line is the total calculated bolometric light curve for case c and the 
arrow marks the time of the BH emergence for this case. 
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where M^h is the BH mass, LEdd = -iTrGMbhc/kes = 1.3 x lO^^(fces/O.4cni2g-i)-i(Affc,i/M0) and 
MEdd = LEdd/c^ = 2.3 X 10~^{kes/0Acm^g^^y^{Mbh/MQ)MQ yr~i are the Eddington luminosity and 
accretion rate and k^s is the electron scattering opacity. 

The late-time evolution of the accretion rate turns out to be in fair agreement with that of the 
hydrodynamic models computed by Colpi, Shapiro & Wasserman (1996; hereafter CSW) for polytropic 
flows with no net radiation flux. In fact, prior to recombination, the flow is radiation dominated and 
very optically thick to emission-absorption so that, from the dynamical point of view, it behaves like an 
adiabatic fluid gas with T = 4/3. On the other hand, by the time that recombination takes place and the 
flow becomes transparent, radiation and pressure forces at around the accretion radius = GMbh/c^ are 
very small and the hydrodynamic evolution turns out to be close to that of a pressureless fluid (dust). 
Introducing the initial accretion timescale tafi = GMbh/c^ o (where Cgfl is the initial gas sound velocity) and 
the expansion timescale to = rout/Vo, CSW have shown that, for envelopes with ia,o/^o < Ij pressure forces 
are important and M evolves along a sequence of quasi-static Bondi-like states, with a slowly decreasing 
density at large distances. However, if the envelope has ta,o/tQ > 1, then the gas behaves like dust, and it 
can be shown that the late-time evolution of M is given by (CSW; equation [29]) 

M 47r2/3 

potokesc — ■ (3) 



Msdd 9 \*o 
Substituting equation (^ into equation (^), we obtain an analytic estimate of the late-time accretion 
luminosity for a pressureless fluid (ZCSW) 

L- 5.7x10 J l,0.4cm2g-ij \Mq 



V 5/6 / , \ 20/9 / , \ -25/18 
_P0 \ / Jo_\ [ J_ 



0.1 gem ^) \lhry \lyr 



erg s ^ . (4) 



We will use equation (||) to describe the evolution of the helium mantle of SN 1997D after the shock 
has passed through the envelope and after the reverse shock (formed at the hydrogen-helium interface) has 
reached the compact remnant (secondary fallback). 



4. Emergence of the black hole in SN 1997D 

The physical conditions of the He mantle after the passage of the shock, the amount of primary fallback 
and the mass of the central BH in SN 1997D depend sensitively on the explosion energy. Successful light 
curve and spectral fitting of SN 1997D show that the energy of the explosion was small (a few x 10^° erg) 
and a significant fraction of the He mantle may have fallen back. Here we assume that, because of significant 
primary fallback, the mass of the BH is M^h — 3A/0 (of which ~ 1.4Mo represents the mass of the iron 
core and ~ 1.6Mq is accreted from the helium mantle). A few hours after the explosion, the internal 
energy in the He mantle of SN 1997D turns out to be typically E ~ 10^^ erg g~^, yielding a sound velocity 
Cs,o = (r£'/3)^/^ ~ 2 X lO"^ cm s^^ (where the polytropic index is F = 5/3 since the flow is gas-pressure 
dominated) and an accretion timescale ta,o ~ 3 x 10^ s. The characteristic value of the expansion timescale 
of the helium mantle inferred from low energy explosion numerical models turns out to be to <^ 1-2 hr (T.R. 
Young; private communication). Since ia,o > to-, the late time hydrodynamic evolution is described within 
the dust approximation. Accordingly, at late times the accretion luminosity of SN 1997D will be given by 
equation (||) with /i = 4/3 and kes = 0.2 cm^ g~^ to account for the helium composition. 
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In Figures [| and |[ we plot the estimated late-time accretion luminosity of the He mantle of SN 1987A 
and SN 1997D. For SN 1997D we plot four curves corresponding to different initial values for the average 
density and the expansion timescale, which span a physically plausible range of the parameter space. In 
order to bracket uncertainties in the value of the average density of the helium mantle a few hours after 
the explosion, we select three values of po corresponding to outer radii in the interval thc — 3 x 10^^-10^^ 
cm. As can be seen from Figure ^, the putative BH in SN 1997D should emerge above the radioactive decay 
emission of the envelope approximately 3 years after the explosion. Since the luminosity by radioactive 
decay falls off exponentially with time, while the accretion luminosity drops only as a power law, this 
estimate is not very sensitive to the actual value of the initial parameters of the He mantle. An estimate of 
the time of the BH emergence above the radioactive decay emission of the envelope can be obtained by 
equating equations (|l|) and (||). We get 



^ In ( ^1 ) = 6.2 - In 



VO.5/ 



-4/3 

X 



0.4cm2g-iy \Mc, 

X 5/6 / ^ 20/9 

Po \ /to 



-1 

X 



0.1 gem J \ 1 hr 



(5) 



Thus, te ~ 6ro ~ 2 yrs, and its value depends only logarithmically on the initial parameters. On the other 
hand, as shown in Figure in the four cases considered, the accretion luminosity at emergence (equation 
Q) varies significantly within the interval 10'^^-5 x 10^^ erg s^^. Thus, uncertainties in the inferred 
parameters of the mantle do not affect significantly the time of the emergence of the accretion luminosity 
above the emission of ^^Co, but they do cause large uncertainties in the magnitude of the luminosity at 
emergence, since L cx pQ^^t^^^^ . Note that, for cases a, b, and c considered in Figure 2, the mass accreted 
during secondary fallback is ^ a few Mq, leading eventually to the formation of a rather massive black hole. 

The possible presence of other radioactive isotopes in the ejecta adds an element of uncertainty to 
the above discussion. As an upper limit in plotting the radioactive decay luminosity of ^''Co and ''^Ti in 
Figure ^ we have adopted the same rate of 7-ray deposition, at an abundance 40 times lower, as in SN 
1987A (~ 1.25 x 1O~^M0; see Section || for a discussion). As can be seen from the Figure, even in the least 
favorable case considered (case d) , the accretion luminosity at emergence is as large as the decay luminosity 
of "Co. 

By comparing Figures 

I and |, it appears that SN 1997D is more promising than SN 1987A for 
detecting an accreting BH because of the lower amount of radioactive elements present in the envelope 
(which shifts the dashed curves downward) and the higher density of the He mantle (which moves the 
dotted curves upward). 



5. Discussion and conclusions 

We have shown that, if a BH has formed in SN 1997D, the late-time accretion luminosity may 
dominate the bolometric light curve ~ 3 years after the explosion. Depending on the actual value of the 
initial density and expansion timescale of the He mantle, the bolometric luminosity at emergence falls in 
the interval — 2 x lO'^^-lO'^'' erg s^^. After breakout, the light curve is expected to fall off as t^^^/^'*. To 
estimate the apparent magnitude at emergence, we assume a black body emission spectrum and extrapolate 
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from the present temperature to get an effective temperature T^ff = 4,000-6,000 K. This is reasonable 
since we expect that the inner accreting region is ionized at T ~ 10^ K and that dust may form in the 
outer ejecta reprocessing the radiation produced by accretion and re-emitting it at a lower temperature. 
Adopting a distance modulus fi = 30.64 mag for NGC 1536, the estimated apparent visual, red and I-band 
magnitudes at emergence are: niy ~ 27, ~ mj ~ 26 if — 10'^^ erg and niv ~ 31, ~ m/ ~ 30 if 
Le — 2 X 10'^^ erg s~^. If Le ^ 2 x 10'^^ erg s~^, the BH could be detectable for a few years in the red and 
I band with the Hubble Space Telescope, until the luminosity drops below threshold. Clearly, the precise 
value of the initial parameters is crucial in establishing the luminosity at emergence and hence whether 
or not the accreting BH can be detected. Even if it is only marginally visible (Le 2 x 10^^ erg s~^), 
subtracting off the ^^Co decay luminosity would reveal a power-law decay accretion luminosity which would 
confirm the presence of a BH. We note that a supernova with properties similar to SN 1997D exploding in 
a nearer galaxy would appear much brighter and could reveal the presence of a newly formed BH from the 
distinct behavior of the late-time light curve. 

The ^-25/18 pQ-y^fgj. la^yj, decay of the light curve is a signature of the presence of a black hole accreting 
from the stellar envelope and follows directly from equations (2) and (3). We note that this result is not 
valid for an accreting neutron star after the accretion shock has reached the trapping radius (Chevalier 
1989). The reasons are that radiation pressure cannot be neglected at so early an epoch (^ few years) and 
the accretion luminosity efficiency is quite different. 

A number of caveats apply to our result. The formation and the actual mass of the BH depend on 
the amount of mass which falls back during the first hours after the explosion. This parameter is quite 
sensitive to the explosion energy, as shown by Woosley & Weaver (1995). In addition, the state of the flow 
in the He mantle is determined by the global hydrodynamics of the explosion. The estimate of the late-time 
accretion luminosity given in Section ^ assumes that the opacity in the accreting region is contributed 
mainly by electron scattering and that the thermal and dynamical effects of radioactive heating can be 
neglected. Early in the evolution, in the inner accreting region, the temperature is relatively high {T ^ 10^ 
K) so that all elements are almost completely ionized. In these conditions, the dominant contribution to 
the Rosseland mean opacity comes from electron scattering and hence our assumption is reasonable. Later 
on, the temperature decreases and bound-free and bound-bound transitions (especially those of heavy 
elements) may be important. A lot of photons could remain trapped in the flow. Analogous considerations 
apply to the expanding envelope at late times. For example, if the mean opacity of the envelope is ~ lO'^fces, 
the time taken by photons to diffuse out of the He mantle would be ~ 20 years for the typical parameters 
of SN 1997D. Finally, in deriving our estimate we have assumed that, after 1-2 years, the amount of heat 
deposited by 7-rays from radioactive decay is not sufflciently large to alter the character of the flow. Given 
the small amount of heavy elements in the ejecta of SN 1997D, this assumption seems reasonable. 

Improved radiation-hydrodynamic calculations are necessary to provide a theoretical guide for 
identifying the nature of the compact remnant. 

It is a pleasure to thank Massimo Turatto, Tim Young and their collaborators for providing us with 
data for the helium mantle of SN 1997D in advance of publication and for several useful discussions. We 
thank the referee for useful comments. This work was supported in part by NSF grant AST 96-18524 and 
NASA grants NAG-5-2925 and NAG-5-3420 at the University of Ilhnois at Urbana-Champaign. 
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